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We revisit the phase transition from the Neel ordered to a valence bond solid (VBS) state in 
the two-dimensional Ji — J2 antiferromagnetic Heisenberg model. In the first part we address the 
question whether or not this transition could be an example of a second order phase transition due 
to a deconfinement of spinons. We give arguments based on series expansion and spin- wave theory 
that this is not the case and the transition is most likely first order. The method proposed here to 
detect first order phase transitions seems to be very sensitive and might be useful in other models 
as well. In the second part we analyze possible VBS patterns in the magnetically disordered phase 
based on numerical data for different susceptibilities, obtained in the ordered phase, which test the 
breaking of lattice symmetries. We conclude that a columnar dimerization pattern is the most likely 
candidate. 
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I. INTRODUCTION 

It is well known that microscopically different systems 
can show similar behavior near a critical point, a phe- 
nomena termed universality. The universal behavior is 
caused by the fact that only a small number of long-wave 
length degrees of freedom are relevant for most physical 
quantities at the critical point. It is indeed often suffi- 
cient to consider an effective theory in which all modes 
other than the order parameter have been eliminated. 
The concept of the order parameter and effective theo- 
ries based on this quantity has been developed by Landau 
and Ginzburgi 

Recently, it has been argued that there is a second 
order phase transition in the S = 1/2 square lattice 
antifcrromagnet between the Neel state and a paramag- 
netic valence bond solid state which is not described by 
a Ginzburg-Landau (GL) type critical theoryi^ Instead 
the proposed theory involves fractional degrees of free- 
dom (spinons) interacting with an emergent gauge field. 
One of the best studied models where this scenario could 
possibly be realized is the spin-1/2 Ji — J2 Heisenberg 
antifcrromagnet in two dimensions 

H=J 1 J2S i S j + J 2 J2 S i S j- (!) 

nn nnn 

Here Ji > is the nearest-neighbor interaction and 
J2 > a frustrating next-nearest neighbor exchange. 
There are two well understood limits: For Ji ^ 0, J2 = 
the model is just the usual two-dimensional Heisen- 
berg antifcrromagnet which is known to possess Neel or- 
der although with an order parameter M « 0.3 which 
is reduced compared to its classical value M = 1/2. 
For J\ — > and J2 / and fixed, on the other 
hand, both sub-lattices are Neel ordered and J\ then 
induces a so called collinear order. For general cou- 
plings Jx,J%, both limiting ground states become frus- 



trated. Therefore a parameter region might be expected 
where the magnetic order vanishes and a spin liquid or 
VBS ground state is formed. Numerical studies includ- 
ing exact diagonalization-, variational Quantum-Monte- 
Carlfj^ as well as series expansionLSiSiiSiii indeed indi- 
cate that for 0.4 < g < 0.6, with g = J2/J1, no magnetic 
order exists. 

To address the question of whether or not the ground 
state in the non-magnetic region breaks a lattice symme- 
try, the response to a field 

F 1 = 5j2(-l) i Si,jSi+i,i (2) 

has been calculated ^ifliiLiSi The series expansion 
studiesi2*ii show that the corresponding susceptibility 
becomes very large or even diverges when g « 0.4 is 
approached from the Neel phase, indicating that transla- 
tional symmetry by one site is broken and a VBS state 
is formed. The results obtained in Ref. by a varia- 
tional quantum Monte Carlo (QMC) technique seem to 
support this scenario. Furthermore, the QMC data have 
been shown to be in good agreement with an exact di- 
agonalization of clusters with TV = 16, 32 spins. In a 
later exact diagonalization study, however, it has been 
shown that the susceptibility decreases when going from 
the 4 x 4 to the 6 x 6-clusteriiS Further evidence in favour 
of a breaking of translational invariance has been ob- 
tained by a dimer series expansion showing directly that 
the corresponding order parameter is indeed nonzero in 
this phased On the basis of the series expansion data 
we therefore believe that the existence of a homogeneous 
spin liquid phase in this parameter region, as proposed 
in Ref. 6, is highly unlikely. 

A very important result of the series calculations in 
Refs. IMlOMlll is that the point g cv , where magnetic order 
vanishes, and the point g c< p, where dimer order becomes 
established, are very close or even coincident. Further- 
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more, a crossing of energies between an Ising expansion 
and a dimer expansion, which would be the indication of 
a first order transition, could not be detected. This lead 
to the assumption that the transition is second order. For 
a second order transition it is, however, difficult to under- 
stand within GL theory why g cv and g cr p should be equal. 
Each phase has a different broken symmetry (spin rota- 
tional symmetry versus lattice symmetry) so that one 
would naively expect that each transition is described 
by its own effective theory containing only the staggered 
magnetization (order parameter in the Neel phase) or 
the dimer order parameter, respectively. In this case the 
transitions should be independent from each other. That 
this seems to be not the case tells us that the two order 
parameters must be related. A GL-type theory should 
therefore at least contain also terms describing the in- 
teraction between the two order parameters. Such an 
effective theory has been proposed in Ref. However, 
it has been found that within this theory g cv and g c( k will 
only be identical if the non-magnetic phase has mass- 
less excitations. This would correspond to a transition 
to a translationally invariant spin liquid phase which we 
believe can be ruled out based on the numerical data 
obtained in Ref. 0. 



The new effective theory for the second order phase 
transiti on b etween the Neel and a VBS state proposed 
in Refs. l2l.'j| follows an entirely different route. Here the 
order parameters in the two phases are represented in 
terms of fractional degrees of freedom (spinons) which 
become deconfined exactly at the critical point. As the 
spinons are the constituents of both order parameters 
this would offer a natural explanation for a direct second 
order phase transition between these at first sight very 
different phases. Among the models proposed to show 
such a transition is the two-dimensional spin- 1/2 model 
with a four-spin exchanged Quite recently, however, it 
has been argued that the transition in this model is more 
likely to be first order^4 In the present work we will 
address the same question about the order of the phase 
transition for the J\ — J2 model. 



Our paper is organized as follows: In section[H]we give 
arguments based on series expansion and spin-wave the- 
ory why the phase transition is most likely a weak first 
order instead of a second order transition. In section 
IIIII we will discuss numerical data for three different sus- 
ceptibilities probing the VBS order in the non-magnetic 
phase. These susceptibilities are obtained in the Neel 
phase where the ground state is known and the series is 
therefore not biased. We will discuss why these data pro- 
vide additional evidence against the deconfinement sce- 
nario and will conclude that the VBS order is most likely 
of the columnar dimer type. The last section presents a 
summary and conclusions. 



II. ORDER OF THE PHASE TRANSITION 

Usually a first order phase transition is detected in 
series calculations by looking for the crossing of energies 
obtained by expansions starting from different states. For 
the J\ — J2 model such an energy-crossing has been de- 
tected between an Ising expansion in the collinear regime 
and a dimer expansion in the VBS phase. - This shows 
that the transition at g ~ 0.62 from the VBS to the 
collinear state is first order (see also Fig. [7J. For the 
transition from the Neel to the VBS state at g ~ 0.4, 
on the other hand, no crossing has been found. More 
precisely, the energies for an Ising expansion and vari- 
ous dimer expansions are so close over a relatively large 
parameter regime around the transition point that it is 
not possible to decide within the accuracy of the series 
if there is a crossing or not (see Ref. and Fig. 0) . In 
all previous series studies it has been implicitly assumed 
that the transition is second order. Here we propose a 
more sensitive method to distinguish between first and 
second order transitions and conclude that the transition 
is most likely weak first order. 

Let us consider the ground-state energy e{5) for the 
Hamiltonian with the field in Eq. J2J included for 
\5\ <C 1 and g < 0.4. We have calculated e(S) for dif- 
ferent g by Ising series expansion. Using an Ising ex- 
pansion means that we start with a state which breaks 
spin-rotational symmetry whereas the lattice symmetries 
are intact. Obviously, it is then impossible to restore 
spin-rotational symmetry and break lattice symmetries - 
as would be required when going from the Neel to the 
VBS state - in any finite order in the expansion. How- 
ever, we can expect that an instability of the state we are 
starting with is signaled by a susceptibility, with respect 
to the corresponding symmetry breaking field, which is 
divergent. 

In Fig. n w e present our numerical data. For a fixed g 




FIG. 1: Ising series data for the ground state energy e(5) with 
the field Fi included and g = 0.0, 0.1, 0.2, 0.3. The lines are a 
guide to the eye. 
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we have fitted e(<5) by a polynomial of the form 

e(5)-e(0) = |5 2 + ^ C - 
The susceptibility is then given by 



G 



Xi 
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As the series for the ground state energy shows better 
convergence than the series for the susceptibility itself 
(which was calculated in Ref . flOf) , we were able to obtain 
Xi with much smaller error bars than before as shown in 
Fig- El The new data are nevertheless consistent with the 
old data within the given error bars. A strong response 




FIG. 2: The susceptibility xi calculated as described in the 
text (solid line) compared to the old data from Ref.^3 (dashed 
line). The dotted line represents a linear fit of the new data 
1/Xi =6.58 -(0.43- g). 

to the field Fi is visible indicating that translational sym- 
metry is broken in the non-magnetic phase. If the phase 
transition with respect to the corresponding order param- 
eter would be second order we expect 1/xi oc (g C — g) 1 * 
where g c( j, is the critical point and 7^ the critical expo- 
nent. In mean-field theory 7^ = 1 and we can indeed 
obtain a nice linear fit of 1/xi as shown in Fig. |2 Note, 
that 70 is not expected to change dramatically even if 
fluctuations are taken into account (an 0(l)-model, for 
example, would have 7^ ss 1.2) so that changes to the 
value of g c $ would be minor. 

To study different possible scenarios for this phase 
transition we consider an effective field theory for the 
magnetically ordered phase. In the effective field the- 
ory for a two-dimensional antifcrromagnet in the ordered 
phase no topological term (Berry phase) is present^ One 
can therefore describe the system by the following 0(3)- 
model 
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H v = ^ + C '(V") 2 + m 2 v v 2 } + ^> 2 ) 2 . (5) 

Assuming a second order phase transition at a critical 
point g cv we have m 2 = a v (g — g C u) 7 " and u v > 0. Here 



7„ is the critical exponent for the staggered magnetic 
susceptibility with j v — 1 in mean-field theory At g < 
g cv the vector field v will then show a nonzer o ground 
state expectation value (v) = \J a v (g cv — g)/u v . 

Consider now the case that we are in the magnetically 
ordered phase and add the field F\ as given in Eq. (0) 
with \8\ <C 1. The Neel order will then coexist with a 
small dimerization described by a scalar field 

If a second order phase transition with respect to <f> at 
a critical point g C( p would occur, we would have = 
o-4>(g c <p — g) 7 * and > 0. We also want to include an 
interaction between the vector and the scalar field. The 
lowest order coupling term allowed by symmetry is 



(7) 



The effective field theory in the ordered phase for 6 ^ 
is then given by H = H v + + H mt and we will have a 
nonzero ground state expectation value 



~ Ar 4 
A 7 



(D(S 7 ) 



(8) 



with A = nifj, + u V( f,(v) 2 . This leads to a ground state 
energy given by 

= o) = -^ s2 ^ s4+ ^^ s6+0 ^ • 

(9) 

This regular expansion in powers of S 2 will exist for 
any g < g cc f, but the radius of convergence will become 
smaller and smaller when the assumed critical point is 
approached. Finally, this expansion will break down and 
directly at the critical point e(<5) will show a scaling with 
a critical exponent which is in general non-integer. Nev- 
ertheless, in the parameter regime where this expansion 
is valid we expect the coefficient of the S -term to be pos- 
itive because > for a second order phase transition. 

We have fitted e(5) shown in Fig. fflby the polynomial 
given in Eq. © but studied so far only the coefficient 
of the quadratic term which gives the susceptibility (see 
Fig. [21 • The coefficient of the quartic term obtained by a 
fit of these data is shown in Fig. Surprisingly, this co- 
efficient is negative. Before we discuss the consequences 
we will check if the coefficient follows the form predicted 
in 0. The coefficient A — l/xi so that we can use the 
linear fit shown in Fig.[3for this parameter. With these 
values for A the data in Fig. [3] are in very good agree- 
ment with Eq. JSJ) where the single parameter fit (shown 
as dashed line) yields = —5.68. 

A check of our series data for b is provided by calcu- 
lating this coefficient in spin-wave theory. We know that 
linear spin-wave theory (LSWT) yields accurate results 
for quantities like ground-state energy or sub-lattice mag- 
netization for the two-dimensional Heisenberg antifer- 
romagnet without frustration^iii Surprisingly, LSWT 
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FIG. 3: The coefficient b of the quartic term in Eq. ob- 
tained from a fit of the data in Fig. (black squares) . The 
dashed line represents a fit b(g) = — 5.68/[6.58(g — 0.43)] 4 . In- 
set: The coefficient b from LSWT and MFSWT in comparison 
to the series data. The lines are a guide to the eye. 



gives also a phase diagram for the J\ — J2 models which 
is very similar to the one found by numerical calcula- 
tions. In particular, a non-magnetic ground state for 
g G [0.38, 0.51] is found. In a spin-wave theory where 
the quartic terms are treated self-consistently in one-loop 
approximation (MFSWT), however, the phase diagram 
changes dramatically and Neel order remains stable up 
to g ~ 0.6^ This tells us that we cannot trust spin- wave 
theory in the strongly frustrated regime. If we consider 
g close to zero, on the other hand, we might expect that 
spin- wave theory results are reasonable. 

In complete analogy to the series calculations we have 
calculated the ground-state energy e(6) in LSWT and 
MFSWT for the J x - J 2 model with the term © included. 
The coefficient of the <5 4 -term obtained from these calcu- 
lations is shown in the inset of Fig. [3] in comparison to 
the series data. Again we find that b is negative and 
decreases with increasing g. Quantitatively, b is larger, 
both in LSWT and MFSWT, as in series and decreases 
more slowly. 

A negative u,), in Eq. © means that the phase transi- 
tion with respect to <ft will be first order and the assumed 
critical point g c $ will never be reached. For 5 = the 
order parameter ((f)} will instead jump from zero to some 
finite value. To determine the transition point one also 
needs to know r^ which has to be positive for stability 
reasons. In principle, r^ can be determined from the 
sixth order term in |||. In practice, the errors in the 
numerical data for e(<5) and in the parameters A^u^ are 
too large to determine reliably. Due to the interac- 
tion term (J7J a jump in (0) will induce a simultaneous 
jump in the Neel order parameter (v). Depending on the 
strength of the interaction u vi p two scenarios are possible: 
The Neel order parameter could jump to zero yielding a 
direct first order transition from the ordered phase to a 
disordered phase with broken translational invariance. A 



more exotic scenario, where (v) jumps to a smaller but 
nonzero value and then decreases further before vanish- 
ing at a critical point, is also possible. This would imply 
that there is a region in the phase diagram where dimer- 
ization and Neel order coexist. Based on the currently 
available numerical data it is impossible to decide which 
scenario is actually realized. We also want to mention 
that our mean-field treatment of the order parameters 
and the use of the expansion © are a posteriori justified. 
The finding of a first order phase transition means that 
we never get to the assumed critical point g c( f, where the 
length scale for fluctuations would diverge and a mean- 
field treatment would therefore be very questionable. It 
also means that the expansion JSJ exists and its radius 
of convergence will be finite everywhere in the ordered 
phase. 



III. VBS ORDER IN THE NON-MAGNETIC 
PHASE 

In this section we provide additional, independent ar- 
guments against the deconfinement scenario and in favor 
of a first order phase transition. We also address the 
question which kind of VBS order is actually realized. 

In Refs. dealing with a possible deconfined criti- 
cal point separating a Neel ordered from a VBS phase 
in a two-dimensional antifcrromagnet it has been implic- 
itly assumed that the VBS order is either of columnar 
dimer type as shown in Fig. 0f a) or of plaquette type as 
shown in Fig. Hfb). The order parameter for columnar 
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FIG. 4: Possible VBS ordering patterns: (a) columnar dimer 
order, (b) plaquette order and (c) plaquette order with every 
second column of plaquettes shifted by one lattice site. 
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dimerization is given by O c = (—iySijSi+i.j whereas 
the order parameter for the plaquette phase can be rep- 
resented as O p = (—l) l SijSi+i t j + (— l)i SijSij+i. As 
discussed in Ref. they can also be interpreted as a single 
complex order parameter where only the phase is differ- 
ent for the two patterns. At a deconfmed critical point, 
however, this phase will only appear as an irrelevant op- 
erator. Therefore both order parameters are expected to 
show power law correlations at such a point and suscep- 
tibilities testing the breaking of lattice symmetries with 
respect to columnar or plaquette order should diverge 
when the deconfmed critical point is approached from 
the magnetically ordered phase. 
First, consider the field 

F2 = S (SijSi+ij — SijSij+i) (10) 

i,3 

which tests if rotational symmetry is broken, which would 
be the case for the columnar dimer state but not for 
the plaquette state. The corresponding susceptibility X2 
shows only a very moderate increase with g and cer- 
tainly no sign of divergence when the possible critical 
point g cc j, ~ 0.43 is approached (see Fig. EJ). We men- 
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FIG. 5: Susceptibility X2 as a function of g calculated by Ising 
series expansion. The line is a guide to the eye. 

tion here that our results for the susceptibilities xi and 
X2 are in qualitative agreement with Ref. y where the 
finite size scaling of these susceptibilities has been inves- 
tigated by exact diagonalization and variational QMC 
methods. Note that it does not come as a surprise that 
there is no quantitative agreement: the variational QMC 
method explores only a small part of the Hilbert space so 
that a convergence to the true value for the susceptibility 
cannot be expected even if arbitrary large clusters could 
be treated. If a good variational wave-function is chosen, 
however, one might expect that the finite size scaling is 
at least qualitatively similar to the true finite size scaling. 
Next, we consider the susceptibility X3 corresponding 



to the field 

which has already been calculated in Ref. ^3- We re- 
produce the result in Fig. HJ1 \3 even decreases when 

30 1 1 
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FIG. 6: Susceptibility X3 as a function of g calculated by Ising 
series expansion. The line is a guide to the eye. 

the non-magnetic phase is approached which seems to 
indicate that the corresponding lattice symmetry, which 
would be broken for the plaquette but not for the colum- 
nar dimer state, is intact. 

The fact that \2 and X3 do not diverge when g — > g c $ w 
0.43 as would be expected for a deconfmed critical point 
provides an argument independent of the considerations 
in section ^ that the phase transition in the J\ — J% 
mod el is not an example for the scenario proposed in 
Refs.EH Putting the arguments given in sectionlTTlaside. 
this alone, however, does not exclude the possibility of 
another kind of second order phase transition from the 
Neel ordered to a state with VBS order. Nevertheless, in 
any second order scenario we should take the fact that xi 
and X3 do not diverge seriously. This means that in this 
case the VBS order cannot be of columnar dimer type 
because then xn should diverge, and not of plaquette 
type because this would lead to a diverging xz- So in 
such a scenario the quest is to find an ordering pattern 
with lattice symmetries which agree with our findings for 
all three susceptibilities. 

A possible pattern is shown in Fig.^c). Translational 
symmetry by one site along the x-axis is broken leading 
to a divergent xi- X3 wm be finite because of the 180° 
rotational symmetry around the axis marked by "1" in 
Fig-EIc) or identical positions, xi Wlu be finite for this 
pattern because of the following symmetry: A rotation 
by 90° around the axis "2" or identical positions, fol- 
lowed by a shift of every second row of plaquettes by one 
lattice site, followed by a shift of every second column 
of plaquettes by one lattice site. Although other pattern 
which have the correct symmetries might be possible, the 
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pattern in Fig.^Jc) seems to be the one with the smallest 
unit cell. 

We have calculated the ground-state energy as well as 
the singlet and triplet dispersion for different g for this 
state by series expansion starting from the decoupled pla- 
quettes shown in Fig. Hfc). The expansion parameter 
x describes the couplings between the plaquettes. Fol- 
ic = the plaquettes are decoupled whereas x = 1 corre- 
sponds to the case we are finally interested in where the 
couplings within and between the plaquettes are of equal 
strength. 

In Fig. [7| we show the ground-state energy obtained 
from this expansion compared to other series data. 




0.3 0.4 0.5 0.6 0.7 

FIG. 7: Ground state energies calculated by different series 
expansions starting from the nonmagnetic states shown in 
Fig. [1] as well as from the magnetic Neel and collinear states. 
The curve 4(a) (open squares) corresponds to the columnar 
dimer state shown in Fig. EI a), the curve 4(b) (black dots) 
to the plaquette state shown in Fig. 21b), and the curve 4(c) 
(red dots) to the plaquette state shown in Fig.^Ic). The black 
triangles show results of Ising series expansions starting from 
the Neel and collinear states. 

The obtained energies for 0.4 < g < 0.5 are similar or 
even slightly lower than the energies obtained from the 
columnar dimer expansion. This shows that the state 
in Fig. QIc) is indeed a possible candidate for the ground 
state in this parameter region. In series calculations it is, 
however, important to make sure that the state which is 
the starting point for the expansion remains stable when 
extrapolating to the isotropic limit x — > 1. We have 
therefore calculated also the singlet and triplet disper- 
sions. The plaquettes for this VBS state are arranged on 



a triangular lattice. With respect to this lattice both sin- 
glet and triplet dispersion show a mininum at momentum 
k = (0,2tt/V3). The series coefficients for the singlet and 
triplet gaps at this point and g = 0.25, 0.45 are given in 
Table [i] Dlog Pade approximants for these series show 
that the singlet gap always vanishes before the triplet 
gap and at a value x c < 1. This indicates that the state 
in Fig. 21c) is unstable. 

For the columnar dimer state, on the other hand, it is 
very difficult to obtain the singlet dispersion by series ex- 
pansion because the state we are starting with consisting 
of decoupled dimers does not contain a singlet excitation. 
We have, however, calculated the triplet gap in this case 
which shows a minimum at k = (0, n) and Dlog Pade 
approximants suggest that in this case the gap remains 
nonzero for x — *■ 1. In addition, the instability of the 
state in Fig. 01c) which manifests itself by a vanishing 
of the singlet gap at k = (0, 2tt/^/3) for x < 1 is an in- 
stability towards a columnar dimerization. Although this 
does not prove that the columnar state is finally stable, it 
makes the columnar dimerization pattern the most likely 
candidate for the VBS order in the non-magnetic region 
of the phase diagram. 



IV. CONCLUSIONS 

We have calculated the ground-state energy e for the 
J\ — J2 model including a small field F\ with strength 
5, which induces a columnar dimerization in the Neel or- 
dered phase, by Ising series expansion. We have argued 
that everywhere except directly at a critical point it is 
possible to expand e(d) in a regular series in 5 2 and that 
this series has a finite radius of convergence which goes 
to zero when the tentative critical point is approached. 
The prefactor of the i5 2 -term in that series gives the sus- 
ceptibility xi w ith respect to iq. The data for this sus- 
ceptibility obtained by an Ising series expansion indicate 
that the Neel state becomes unstable for g > g c 0.43 
and that the ground state for g > g c breaks translational 
symmetry by one site and therefore seems to be of the 
VBS type and not a spin liquid. Based on a mean-field 
treatment of an effective field theory describing the Neel 
state coexisting with the small dimerization induced by 
iq we have argued that the sign of the <5 4 -term in the 
expansion of e(S) determines whether the transition with 
respect to the VBS order parameter is first or second or- 
der. We believe that this is an in general more sensitive 
and less biased method to distinguish between a first and 
a second order transition than looking for a crossing of 
energies obtained by different expansions. From the se- 
ries data we found that the <5 4 -term has a negative sign 
and we showed that the same is true in spin-wave theory. 
Within the presented GL-type theory this means that the 
transition is expected to be first order. Our mean-field 
treatment of the order parameter is a posteriori justified 
because a critical point where such a treatment would 
break down is never reached. 
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TABLE I: Series coefficients x n for the minimum singlet gap A s and triplet gap At for the plaquette state from Fig. 01c) and 
g = 0.25 (g = 0.45), respectively. 



n A s /Ji (g = 0.25) A t /Ji (g = 0.25) A 3 /Ji (g = 0.45) At/Ji (g = 0.45) 

1.500000000 1.000000000 1.100000000 1.000000000 

1 0.000000000 -8.33333334X10" 1 0.000000000 -4.33333334x 10" 1 

2 -4.356195887xl0 _1 -5.722808442xl0 _1 -4.151003339x 10" 1 -5.005966374xl0 _1 

3 -3.880918282 xlO -1 1.4032850965 x 10 _1 -1.832689336x 10 _1 8.2113320127x 10~ 2 

4 -4.828163654xl0 -1 -2.265506840xl0 _1 -1. 537390367x10^ -3.080455443 x 10 -2 

5 -6.032418702 xl0 _1 7.5790435073 x 10" 3 -1.592497102xl0 _1 -1.108096117x 10 _1 

6 -7.988310579xl0 _1 -4.517676083xl0" 2 -1.653261118xl0 _1 1.6084176410x 10~ 2 



In the second part we gave arguments in favor of a first 
order transition which are independent of any effective 
field theory by analyzing two additional susceptibilities 
testing different lattice symmetries. These susceptibili- 
ties were calculated based on an Ising expansion in the 
Neel phase so that the series is not biased by any assumed 
dimerization pattern. We argued that at a deconfincd 
critical point all three susceptibilities considered here are 
expected to diverge and that the fact that \i an d X3 do 
not diverge excludes this scenario. We further argued 
that in any second order scenario the non-divergence of 
%2 would mean that the VBS state is not of the colum- 
nar dimer type and the non-divergence of X3 would mean 
that the VBS state is not of plaquette type either. For 
an assumed second order transition we have been able to 
find a VBS pattern which does have the correct lattice 
symmetries to explain our data for all three susceptibil- 
ities. Series expansion data starting from this pattern, 
however, have proven that this state is unstable with re- 
spect to the columnar dimerization pattern. 

Taking the arguments given in the two parts together 



shows that the transition from the Neel state to a VBS 
state with columnar dimerization is most likely first or- 
der. For a first order phase transition we still expect 
that the susceptibility xi when calculated by Ising se- 
ries expansion diverges because the field F\ directly tests 
the instability of the Neel state with respect to columnar 
dimer order irrespective of the order of the phase tran- 
sition. For the susceptibility Xii on the other hand, we 
can only expect that it diverges if the transition is sec- 
ond order. The fact that \2 does increase with g without 
diverging at g ~ 0.4 therefore further supports that the 
transition is most likely first order and the VBS state 
most likely of the columnar dimer type. 
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